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CHAPTER 1 
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In these lectures I describe long scale properties of fluctuating polymerized mem- 
branes in the presence of network anisotropy and random heterogeneities. Amaz- 
ingly, even infinitesimal amount of these seemingly innocuous but physically im- 
portant ingredients in the membrane's internal structure leads to a wealth of 
striking qualitatively new phenomena. Anisotropy leads to a "tubule" phase that 
is intermediate in its properties and location on the phase diagram between pre- 
viously discussed "crumpled" and "flat" phases. At low temperature, network 
heterogeneity generates conformationally glassy phases, with membrane normals 
exhibiting glass order analogous to spin-glasses. The common thread to these 
distinct membrane phases is that they exhibit universal anomalous elasticity, 
(singularly length-scale dependent elastic moduli, universal Poisson ratio, etc.), 
driven by thermal fluctuations and/or disorder and controlled by a nontrivial 
low-temperature fixed point. 

1. Preamble 

The nature of a membrane's in-plane order, with three well-studied universality 
classes, the isotropic hquiflhexaticliquicElandsolicP, crucially affects its confor- 
mational properties. The most striking illustration of this (discussed in lectures by 
Yacov Kantor and by David Nelson) is the stabilization in polymerized membranes 
(but not fluid ones, that are always crumplecP beyond a persistence lengthPJ) of a 
"flat" phasd^, with long-range orientational order in the local membrane normald^El 
that is favored at low temperature over the entropically preferred high-temperature 
crumpled state. In a beautiful "order from disorder" phenomenon, a subtle interplay 
of wild thermal fluctuations with nonlinear membrane elasticity (made possible by 
a finite shear modulus of a solid) infinitely enhances membrane's bending rigid- 
ity, thereby stabilizing the flat phase against these very fluctuations. A univer- 
sal fluctuation-driven "anomalous elasticity" characterizes the resulting flat phase, 
with length-scale dependent elastic moduli, non-Hookean stress-strain relation, and 
a universal negative Poisson ratio! 1 11 

Given such qualitative distinction between liquid and solid membranes, it is per- 
haps not too surprising that other in-plane orders can have important qualitative 
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effects on membrane's long-scafe properties. In these lectures I will discuss two such 
ingredients, namely, membrane network in-plane anisotropy (Section 2) and hetero- 
geneity (Section 3) in fluctuating polymerized membranes, and will show that these 
seemingly innocuous generalizations, indeed lead to a wealth of new phenomena 
that are the subject of these lectures. 

2. Anisotropic Polymerized Membranes 
2.1. Motivation and introduction 

In addition to the basic theoretical motivation, the interest in anisotropic polymer- 
ized membranes is naturally driven by a number of possible experimental realiza- 
tions, some of which are illustrated in Fig^ One example is a tethered membrane 



Fig. 1. Examples of anisotropic polymerized membranes with (a) a polymerized-in lipid tilt order 
in a bilayer membrane, (b) a weakly crosslinked array of linear polymers, and (c) a spontaneous 
in-plane ncmatic order in a nematic elastomer membrane. 

made through photo-polymerization of a fluid phospholipid membrane exhibiting 
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lipid tilt order. On average the lipids are tilted relative to the membrane normal, 
inducing a vector in-plane order and an intrinsic elastic anisotropy that can in prin- 
ciple be aligned with an electric or a magnetic field and polymerized in. Another 
possible method 1 ^ of fabricating polymerized sheets is by cross-linking a stretched 
out, aligned array of linear polymers, that would clearly lead to an intrinsically 
anisotropic tethered membrane. One other promising candidate is a two-dimensional 
sheet of a nematic elastomeiP] a material that received considerable attention re- 
cently because of its novel elastic and electro-optic properties.El These cross-linked 
liquid-crystal polymer gels exhibit all standard liquid-crystal phases and therefore 
in their nematic state are highly anisotrop ic. 

Almost 10 years ago, it was discovered^! that in-plane anisotropy has a dra- 
matic qualitative effect on the global phase diagram of polymerized membranes. 
As illustrated in Fig|2 it leads to an entirely new "tubule" phase of a polymerized 
membrane (see Fig|3Jl, that is crumpled along one and extended along the other 
of the t wo m embrane axes, with wild undulations about its average cylindrical 



geometry? 



_|_- TUBULE 

4>j>0, <b y =0 



CRUMPLED 






Fig. 2. Phase diagram for anisotropic tethered membranes showing crumpled, tubule and flat 
phases as a function of reduced temperatures t± oc T — , t y oc T — . 



This is possible because in an anisotropic membrane, where symmetry between 
x± and y axes is broken, (e.g., curvature moduli k x± ^ n y ) the x±- and y-directed 
tangents £j_ = d±r and C, y = d y f will generically order at two distinct temperatures 
~ k x± and ~ n y> respectively, thereby allowing two intermediate tubule 
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states: (1) a^-extended tubule with (Cl) ^ 0, (Q = 0, for T% < T < and (2) 
y-extended tubule with (Cl) = 0, (Q ^ 0, for <T < T». 

The direct crumpling transition occurs in such more generic anisotropic mem- 
brane only for that special set of cuts through the phase diagram (like P2) that 
pass through the origin and is in fact tetra-critical. Generic paths (like Pi) will ex- 
perience two phase transitions, crumpled-to-tubule, and tubule-to-flat, that are in 
distinct universality classes. The tubule phase is therefore not only gcncrically pos- 
sible, but actu ally unavoidable, in membranes with any type or amount of intrinsic 
anisotropy\^ 

As illustrated in Fig|2|the tubule is characterized by its thickness Rq, (the radius 
of gyration of its crumpled boundary), and its undulations h rms transverse to its 
average axis of orientation. Quite generally, (for a y-extended tubule of an L±_ x L y 
membrane) they obey the scaling laws: 



Rg(L±, L y ) = L±S R\L y I XjJ 



,(L ± ,L y ) = L<S h (L y /L* ± ) 



(1) 



where the roughness exponent C = v / z i and the anisotropy exponent z = (1 + 
2i/)/(3 — n K ) < 1 are expressed in terms two independent universal exponents: the 
radius of gyration exponent v < 1 and the anomalous elasticity exponent n K for the 
tubule bending rigidity defined by k 



L^ K . The universal scaling functions S^h (x) 



rms 




Fig. 3. Schematic of the y-tubulc phase of anisotropic polymerized membrane, with thickness Rq 
and roughness h rms . 



have the limiting forms: 



Sr{x) oc 



X-v p /z 



for 



constant, for . 



(2) 



S h {x) oc 



constant, for : 
for x 



12 



(3) 



where v p is the radius of gyration exponent of a coiled linear polymer « 3/5. These 
asymptotic forms emerge fr om g eneral requirements (supported by detailed renor- 
malization group calculations ) that in the limit of a very thin tubule (in the curled 
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up _L direction) the tubule's height undulations h rms (L y ) reproduce statistics of a 
linear polymer of length L y and width L±. And, that in the opposite limit of van- 
ishing L y , scaling functions must recover the radius of gyration Rq(L±) ~ L v f of a 
linear polymer of length L± and thickness L y . 

The scaling forms, Eq[3 and 03 imply that for a " roughly square" membrane - 
that is, one with L± ~ L y = L - in the limit L — > oo 

R G (L ± ~ L y = L) oc V , h rms (L ± ~ L y = L) cx L 1 -^ 2 . (4) 

Combining this with detailed renormalization group calculations, that find a strictly 
positive r] K shows that h rms « L for a large roughly square membrane. Thus, 
the end-to-end orientational fluctuations 9 ~ h rms /L oc L^^^l 2 — ► as L — ► oo 
for such a roughly square membrane, proving that tubule order (which requires 
orientational persistence in the extended direction) is stable against undulations of 
the tubule embedded in d = 3 dimensions. 

The rest of Sectionals devoted to developing a model of an anisotropic polymer- 
ized membrane and using it to study phase transitions into and out of the tubule 
phase and the anomalous elastic properties of a fluctuating tubule summarized 
above. 



2.2. Model 



A model for anisotropic membranes is a generalization of the isotropic model dis- 
cussed in David Nelson's lectures.® As there, we characterize the configuration of 
the membrane by the position r(x) in the d-dimensional embedding space of the 
point in the membrane labeled by a D-dimcnsional internal co-ordinate x, with 
d = 3 and D = 2 corresponding to the physical case of interest. 

Rotational and translational symmetries require that the Landau- Ginzburg free 
energy F is an expansion in powers of tangent vectors £± tV = d± iy f and their 
gradients 



F[f(x 

+ t± {d^Y + t y (cyf -f 



K -L (d±r) + K y {dyr) + nj_ y d y r- d\r 



+ - I d D x I d D x' 6^ (f(x) - r(x')) 



(5) 



where we have taken the membrane to be isotropic in the D — 1 membrane di- 
rections Xj_, orthogonal to one special direction, y. While at first glance F might 
appear quite formidable, in fact (aside from the last term), in terms of the tangent 
order parameter £_L, y it has a standard form of the Landau's </> 4 theory. The first 
three terms in F (the k terms) represent the anisotropic bending energy of the mem- 
brane. The elastic constants t± and t y are the most strongly temperature dependent 
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parameters in the model, changing sign from large, positive values at high temper- 
atures to negative values at low temperatures. The u and v quartic elastic terms 
are needed to stabilize the membrane when one or both of the elastic constants t± , 
t y become negative. The final, b term in EqEl represents the self-avoidance of the 
membranes, i.e., its steric or excluded volume interaction. 



2.3. Mean-field theory 

In mean-field theory, we seek a configuration r(x) that minimizes the free energy, 
EqEl The curvature energies k± (p\r) 2 and k v {dyr) 2 are clearly minimized when 
r(x) is linear in x 

r(x) = (Cxx_ L ,C„J/,0,0,...,0) . (6) 

Inserting EqEJinto Eq|SJ and for now neglecting the self-avoidance, we obtain the 
mean-field free energy density for anisotropic membranes 

/mft = \ [tyC 2 v +t±(D- l)Ci + \u' xl _{D - l) 2 Cl + \u m $ + v ±y (D - i)cic, 2 ] , 

(7) 

Minimizing the free energy over order parameters £_l and Cy yields two possibl e ph ase 
diagram topologies, depending on whether u'j_ ± u yy > Vj_ y or u' ±J _u vy < U** 

For u' ± j_u yy > Vj_ y , the phase diagram is given in Fig|2 For t± > and t y > 0, 
average tangent vectors Cj_ and ( y both vanish, describing a crumpled (collapsed to 
the origin) membrane. 

In the regime between the positive i^-axis and the t y < part of the 
t y = {u yy /v±y)t±_ line, lies the y-tubule phase, characterized by £_l = and 
Cy = y/\ty\/ u yy > ; the membrane is extended in the y-direction but crumpled 
in all D — 1 _L-directions. 

The _L-tubule phase is the analogous phase with the y and _L directions reversed, 
C y — and = \/\t±\/uj_± > 0, and lies between the t±_ < segment of the line 
t y = (vj_ y /u'j_j_)t± and the positive t y axis. 

Finally, the flat phase, characterized by both 

1 /2 

Cl = [{\ti_\uyy - \ty\v ±y )/(u' ±± u yy - v ±y )] > , (8) 
Cy = [(\ty\ux± - \t±\v± y )/{u' ±± u yy - vly)] 1/2 > , (9) 

lies between the t±_ < segment of the line t y — (u yy /v± y )t± and the t y < 
segment of the line t y — (v± y /u' ±± )t±. 

For u'^^Uyy < v^_ y , the flat phase disappears, and is replaced by a direct first- 
order transition from _L-tubule to y-tubule along the locus t y — (v± y /u' ±± )t± (sec 
Fig0J . The boundaries between the tubule and the crumpled phases remain the 
positive t y and t± axes, as for u'j_ ± u yy > v\ y case. 
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Fig. 4. Phase diagram for anisotropic tethered membranes showing a new tubule phase, for 
the range of elastic parameters when the intermediate flat phase disappears. A first-order phase 
transition separates y- and ±-tubule phases. 



2.4. Fluctuations and self-avoidance in the crumpled and flat 
phases 

Although anisotropy has dramatic effects on the phase diagram of poly meri zed 
membranes, it does not modify the nature of the crumpled and flat phases!^ To 
see this for the crumpled phase note that when both t± and t y are positive, all 
of the other local terms in Eq[Sl i.e., the k, u, and v-terms, are subdominant at 
long wavelengths. When these irrelevant terms are neglected, a simple change of 
variables xj_ = x'Wt±Jty makes the remaining energy isotropic. Thus, the entire 
crumpled phase is identical in its scaling properties to that of isotropic membranes. 
In particular, the membrane in this phase has a radius of gyration Rg(L) which 
scales with membrane linear dimension L like L v , with v = (D + 2) / (d ',+ 2) in Flory 
theory, and very similar values predicted by e-expansion techniques 

Fluctuations in the flat phase can be incorporated by considering small devia- 
tions from the mean-field conformation in EqEl 



r(x) = [(_LX± + u_L(x),( y y + Uy(x),h(x) 



(10) 



where u a (x) are D in-plane phonon fields and hi (x) are d c = d — D out-of-planc 
height undulation fields. Inserting this into free energy, Eq[5l with t± and t y both in 
the range in which the flat phase is stable, we obtain the uniaxial elastic energy stud- 
ied by TonerPS As he showed, amazingly, fluctuation renormalize this anisotropic 
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elastic energy into the isotropic membrane elastic energy studied previouslyP^EEl 
Therefore, in the flat phase, and at sufficiently long scales, the anisotropic mem- 
branes behave exactly like isotropic ones. 

2.4.1. Anomalous elasticity of the flat phase 

This in particular implies that the flat phase of anisotropic membranes is stable 
against thermal fluctuations even though it breaks continuous (rotational) symme- 
try and is two-dimensional!^ As in isotropic membranes, this is due to the fact that 
at long wavelengths these very t hermal fluctuations drive the effective (renormal- 
ized) bend modulus k to infinitjEEEEl thereby suppressing effects of these same 
fluctuations that seek to destabilize the flat phase, resulting in an Esher-like "order- 
from-disorder" phenomenon. 

Specifically, n(q) becomes wavevector dependent, and dive rge s like q~ llK as q — > 
0. In the flat phase the standard Lame coefficients fi and A are also infinitely 
renormalized and become wavevector dependent, vanishing in the q — > limit as 
fi(q) ~ X(q) ~ q r,u . The flat phase is furt hermore novel in that it is characterized 
by a universal negative Poisson raticP^ which for D = 2 is defined as the long 
wavelength limit q — > of a = \(q) / (2[i(q) + X(q)). The transverse undulations in 
the flat phase, i.e. the membrane roughness h rms scales with the internal size of 
the membrane as h rms ~ L^, with £ = (4 — D — r] K )/2, exactly. Furthermore, an 
underlying rotational invariance imposes an exact Ward identity between rj K and 
i]u, f]u + 2?7 K = 4 — D. This leaves only a single independent exponent, characterizing 
the properties of the flat phase of even anisotropic membranes. 

To appreciate how exotic and unusual this anomalous elasticity really is one 
only needs to observe that most ordered and disor dere d states of matter (systems 
with quenched disorder being prominent exceptions 1 ^), are in a sense trivial, with 
fluct uat ion s ab out them described by a harmonic theory controlled by a Gaussian 
fixed pointP^That is, generically, qualitatively important effects of fluctuations are 
confined to the vicinity of isolated critical points, where a system is tuned to be 
"soft", and characterized by low energy excitations. However, there exists a novel 
class of systems, with flat ph ase of polymerized membranes as a pr ominent mem- 
ber (that also includes smecticP^Eni an d columnar liquid crystalJ ^t vortex lattices 



in putative magnetic superconductors'^, and nematic elastomer Jl^ | 34 | 35 | 36 | 37 | 
whose ordered states are a striking exception to this rule. A unifying feature of these 
phases is their underlying, spontaneously broken rotational invariance, that strictly 
enforces a particular "softness" of the corresponding Goldstone mode Hamiltonian. 
As a consequence, the usually small nonlinear elastic terms are in fact compara- 
ble to harmonic ones, and therefore must be taken into account. Similarly to their 
effects near continuous phase transitions (where they induce universal power-law 
anomalies), but extending throughout an ordered phase, fluctuations drive nonlin- 
earities to qualitatively modify such soft states. The resulting strongly interacting 
ordered states at long scales exhibit rich phenomenology such as a universal nonlo- 
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cal elasticity, a strictly nonlinear response to an arbitrarily weak perturbation and 
a universal ratio of wavevector-dependent singular elastic moduli, all controlled by 
a nontrivial infrared stable fixed point illustrated in FigEl 



anomalous phase — ► 



T* T c 



Fig. 5. Renormalization group flow for anomalously elastic solids, with T c the transition temper- 
ature to the ordered state and T„ a nontrivial infrared stable fixed point controlling properties of 
the strongly interacting ordered, critical-like state. 



2.4.2. SCSA of the flat phase 

The study of such anomalous elasticity in polymerized membranes was initiated by 
Nelson and Peliti using a simple one-loop self-consistent theory, that assumed a non- 
renormalization of the in-plane elastic Lame' constants.^ They found that phonon- 
mediated interaction between capillary waves leads to a divergent bending rigidity 
with rj K = 1 and membrane roughness exponent £ = 1/2. Subsequent controlled 
e = 4 — D and 1 / d c expansions confirmed existence of anomalous elasticity, 
and discovered an additional effect of softening (screening) of in-plane elasticity by 
out-of-plane undulations, that lead to a vanishing of long-scale Lame' coefficients 

Kq) ~ Kq) ~ <P n - 

Here we treat anomalous elasticity within the so-called self-consistent screening 
approximation (SCSA), first applied to the study of the flat phase of polymerized 
membranes by Le Doussal and Radzihovsky®. The attractive feature of SCSA is 
that it is becomes exact in three complementary limits. By construction, it is exact 
in the large embedding dimension d c — > oo limit and agrees with the systematic 
1/d analysis of Guitter, et al0 Because of Ward identities associated with rota- 
tional invariance it is exact (at any d) to lowest order e — 4 — D, i.e., agrees with 

[71 

one-loop results of Aronovitz and Lubensky. (It would be very interesting to check 
predictions of SCSA to order e 2 , however technical difficulties with two- loop cal- 
culations have so far precluded thiJS^l). And finally it gives exact value of r) K (D) 
for d c = 0. Given these exact constraints, it is not surprising that for physical 
dimensions (D = 2, d = 3) SCSA exponents and the universal negative Poisson 
ratio® compare so well w ith latest, largest scale numerical simulations, discussed in 
lectures by Mark Bowick.^0 

As discussed above, at sufficiently long scales, flat phase of an anisotropic poly- 
merized membrane is identical to that of an isotropic onJ23 and is described by 
an isotropic effective free energy that is a sum of a bending and an in-plane elastic 
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energies: 

F[h,u] = J d D x[^(V 2 h) 2 + f m 2 a/3 + ~ U 2 aa } , (11) 

where the strain tensor is u a p = ^(d a U0 + dpu a + d a h ■ dph). To implement the 
SCSA in the flat phase it is convenient to first integrate out the noncritical phonons 
fields u a obtaining a quartic theory for interacting Goldstone tangent vector modes 
d a h, described by free energy 

F[h] = J d D x[^{V 2 h) 2 + ^(dj ■ d h)R aPnS (d-yh ■ d s h)} (12) 

where for convenience, we rescaled Lame coefficients so that the quartic coupling is 
of order l/d c . The four-point coupling fourth-rank tensor is given by 



where P?g = 8 a p — q a qp/q 2 is a transverse (to q) projection operator. The con- 
venience of this decomposition is that K = 2/i(2/i + DX) / (2/i + A) and fi moduli 
renormalize independently and multiplicatively. _____ 
To determine the renormalized elasticity, we set up a l/d c -expansiorP --- 
for 2-point and 4-point correlation functions of h, and turn them into a closed self- 
consistent set of two coupled integral equations for self-energy E(k) that define 
SCSA 

E(k) = ^-kakpk^ks J -R Q/3i7(5 (q)G(k - q) , (14) 
Rup^siq) = i? Q /3, 7 5(q) - i?a/3 :M iy(q)n AIlyiM / I ,/(q)i? M / l/ / i7( 5(q) , (15) 



where G(k) = l/(K,fc 4 + S(k)) = l/(«;(k)fc 4 ) is the renormalized propaga- 
tor, i? Q/ 3, 7 «5(q) the bare quartic interaction vertex and -R a /3,7<5(q) the "screened" 
interaction vertex dressed by the vacuum polarization bubbles Tl a p^s(c() = 
J p PaP/3P~tP8G(p)G(ci — p), and tensor multiplication is defined above. 

In the long wavelength limit these integral equations are solved exactly by a 
renormalized propagator G(k) w 1/E(k) w Z/fc 4- ''", with Z a non-universal am- 
plitude. Substituting this ansatz into Eas ll4ll~ and solving for the renormalized 
elastic Lame moduli, we find that indeed they must vanish as a universal power, 
with fi(q) cx A(q) ~ q nu , and the phonon anomalous exponent r\ u = 4 — D — 2i] K 
related to tTk by dimensional analysis (power-counting on q) first obtained by Nelson 
and PelitiP 'This recovers the celebrated exponent relation enforced by the underly- 
ing rotational invariance of the membrane in the embedding space. The information 
about remaining independent exponent r] K resides in the ^-dependent amplitudes 
of above equations. Cancelling out the nonuniversal scale Z, we obtain 

ii« 'r[iD + i,Jr[2-, K -iD]r[,,, + iD]r[iD + 2-i,J' 
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that determines r] K {D^d) within SCSA. For D = 2 this equation can be simplified, 
and one finds 

"- (D = 2 ' 4) = Wl6-24 + ^ <17) 
Thus for physical membranes (d c = 1) we obtain: rj K = 0.821, rj u = 0.358 and: 

C = 1 - — = 1 = 0.590 . (18) 

2 V15 + 1 

At long length scales SCSA also gives a universal ratio between renormalized in- 
plane elastic moduli determined by Rap^/sio), Ea ll5l and therefore predicts a uni- 
versal and negative Poisson ratio 

a ^ 2M g )+(g-DA( g ) = 4- (l9) 

that compares extremely well with the most recent and largest simulations!^ 
Expanding the result EqEl m l/d c we obtain: 

8 D - 1 Y[D] 1 

7?K_ 4£> + 2r[fpr[2-f] + ° { dl h (20) 

= 1 + (21) 

which coincides with the exact resulJ^EI as expected by construction of the SCSA. 
Similarly, expanding Ea llBl to first order in e = 4 — D we find: 

* = ^ri?^ ' (22) 

also in agreement with the exact result.^ This is not, however, a general property 
of SCSA, but is special to membranes, and can be traced to the convergence of the 
vertex and box diagrams. 

Because in the flat phase, widely intrinsically separated parts of the membranes 
(i.e., points x and x', with |x — x'| large) do not bump into each other (i.e., never 
have r(x) = r(x')), the self-avoidance interaction in Eq0is irrelevant in the flat 
phase. Hence we expect above predictions to accurately describe conformational 
properties of physical isotropic and anisotropic polymerized membranes. 



2.5. Fluctuations in "phantom" tubules 

Thermal fluctuations of the (y-) tubule about the mean- field state r G (x) = Cy{v-> 0) 
are described by conformation, 

r*(x) = {C v y + u(x), £(x)) , (23) 

where h(x) is a d — 1-component vector field orthogonal to the tubule and it(x) is 
a scalar phonon field along the tubule (taken along y-axis). The order parameter 
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C, y is a tubule extension scale that is slightly modified by thermal fluctuations from 
the mean- field value given in Sec 12. 31 and is determined by the condition that all 
linear terms in /i(x) and u(x) in the renormalized elastic free energy exactly vanish. 
This criterion guarantees that h(x) and u(x) represent fluctuations around the true 
tubule ground state. 

Inserting the decomposition Eq[2Blinto the free energy EqO neglecting irrelevant 
terms (e.g., phonon nonlinearities), and, for the moment ignoring the self-avoidance 
interaction, gives F = F m f t + F e i , where F m [ t is given by Eq[7| (with £j_ — and 
F e i [u(x), h(pc)] is the fluctuating elastic free energy part 



Fa = \ I <l D : ■'■L'l.'i 



K{d 2 y h) 2 + t(d^h) 2 + g±{d^uf + g y (d y u + l -{d y h) 2 ) 2 



, (24) 



where k = K y , t = t± + v ±y C 2 , g y = u yy ( 2 /2, g ± = t+ Uj_ y Cy, and 7 = t y + u yy Q 
are elastic constants. 

The underlying rotational invariance of the tubule in the embedding space is 
responsible for two essential features of F e i. Firstly, it enforces a strict vanishing 
of the (d y h) 2 tension-like term, with curvature (d 2 h) as the lowest order harmonic 
term in the Goldstone tangent mode d y h. The result is highly anisotropic bulk 
elastic energy. Rotational invariance also ensures that nonlinear elasticity can only 
come in through powers of nonlinear strain tensor E(u,ti) = d y u + ^(d y h) 2 , and 
that this property must be preserved upon renormalization. 

2.5.1. Anomalous elasticity of the tubule phase 

Within harmonic approximation tubule, bulk rms transverse height undulations are 
given by 

"~ 1 q±dq y 1 _ t 5/2-d 

(25) 



This suggests that for "phantom" tubules, the u pper critical dimension D uc — 
5/2, contrasting with D uc = 4 for the flat phaseMJ^I This also implies that rms 
fluctuations of the tubule normals are given by 

(\Sn y (x)\ 2 )cxL 3 / 2 - D . (26) 

Since this diverges in the infra-red Lj_ — > 00 for D < Di c = 3/2, this harmonic 
bulk mode analysis (ignoring anomalous elasticity and zero modes) suggests that 
the lower critical dimension Di c below which the tubule is necessarily crumpled is 
given by D ic = 3/2. 

As for the flat phase, one needs to assess the role of elastic nonlinearities that 
appear in free energy F e i , Eq|^] One way to do this is to integrate out the phonon 
u (which, at long scales can be done exactly). This leads to the only remaining 
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nonlinearity in h 

F**h[h] = \ I d D - x xdy(d y h ■ d y h)V h (d y h ■ dyh)}, (27) 
with the Fourier transform of the kernel given by 

VhW = 9 ff l a • (28) 

Because of the |q|j_ « q 2 (z — 1/2) anisotropy of the bulk h modes, V^(q) scales 
as g±q± at long scales, and therefore is strongly irrelevant near the Gaussian fixed 
point as long as g± is not renormalized (but see below). It is straightforward to 
verify to all orders in pert urb ation theory, that in a phantom tubule, there is no 
such renormalization of <7j_P^ 

However, as asserted earlier, the full elasticity EqE] before u is integrated out, 
is anomalous, because g y (q) is driven to zero as q y — » 0, according to 

9 v (cL)=q^S g (qy/ql), (29) 

with S g (x) universal scaling function: 

, . / constant, x -> oo 

and exact exponents: 

z = i , r)u = 5 - 2D , (31) 

One simple way to see this is to note that rotational invariance enforces graphical 
corrections to preserve the form of the nonlinear strain tensor E(u, h) = d v u + 
Tj(d y h) 2 . This leads to a relation between r) K , r\ u and the anisotropy exponent z 

2r ]K + r lu = 4-(D-l)/z (32) 

which, together with the defining relation z = 2/(4 — r] K ) reduces to 

2 Vu - (D - 5)r?« = 10 - 4D , (33) 

and for r\ K — gives z and r\ u in Eg 1311 This is supported by a detailed self- consistent 
one-loop perturbative calculation of g y (q) and by direct RG analysis^! 
For phantom membranes with D = 2, rj u = 1 and z — 1/2, we find: 



I y/Ql, % « ^/qT ■ 
This leads to phonon rms fluctuations given by 

( W (x) 2 ) =L 1 / i S u (L y /L 3 / i ) , (35) 
with universal scaling function having limiting form 

/ constant, x -> oo 
S u (x)<xi x/a x ^ Q ^ (36) 
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For roughly square membranes, L y ~ Lj_ = L, so, as L — > oo, L y /L 3 / 4 — > oo this 
gives 

Mx) 2 ) - LT , (37) 
a result that is consistent with simulations by Bowick, et al.. 

2.5.2. Zero-modes and tubule shape correlation 

The tubule's shape is characterized by 

R 2 G = (\h(L ± ,y)-h(0 ± ,y)\ 2 ) , h 2 rms = (\h(x x , L v ) - h(x x , 0)| 2 ) . (38) 

As illustrated in Fig|31 Rg measures the radius of a typical cross-section of the 
tubule perpendicular to its extended (y-) axis, and h rms measures tubule end-to- 
end transverse fluctuations. 

For a phantom tubule these are easily computed exactly using tubule free energy, 
Ea 1241 One important subtlety is that one needs to take into account "zero modes" 
(Fourier modes with qj_ or q y = 0), that, because of anisotropic scaling q± ~ q 2 of 
the bulk modes can dominate tubule shape fluctuations. 

For Rq one finds 

~k B T f d D - l q±_ 1 



R l G = 2{d- D) 



k B T 



L y J L -i (2^-i tq\ 

d D - 1 q ± dq v (1 - e^ ^) 



(^)° tq\ + K (q)g 4 



(l_ e ^-^) 

(39) 



with the first and second terms coming from the q y — "zero mode" and the 
standard bulk contributions, respectively. From its definition, it is clear that the 
= "zero mode" does not contribute to Rq- Standard asymptotic analysis of 
above integrals gives 

R G (L ± ,L y ) = LlS R (L y /L z x ) (40) 
with, for phantom membranes, 

5 - 2D 1 



4 2 ' 

and the limiting form of the universal scaling function Sr(x) given by 



(41) 



S (x) oc I l ^ for X ~* ° f42) 

\ constant, for x — > oo . 

This gives 

R G oc L\ oc L 1/4 , for D = 2 , (43) 

for the physically relevant case of a square membrane Lj_ ~ L y ~ L — > oo, for 
which L y >> Lj_, with bulk mode contribution dominating over the q y = zero 
mode. This result is in excellent quantitative agreement with simulations of Bowick, 
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et all^l who found v = 0.24 ± 0.02, in D = 2. It would be interesting to test the 
full anisotropic scaling prediction of Eas l4QI42l by varying the aspect ratio of the 
membrane in such simulations. For instance, for fixed L± and increasing L y these 
predict no change in Rq- The same should hold if L y is decreased at fixed L±: Rq 
should remain unchanged until L y ~ s/L_l, at which point the tubule should begin 
to get thinner (i.e. Rq should decrease). 

Equations 1401 and 1421 also correctly recover the limit of L y — constant << L\ — > 
oo, where the tubule simply becomes a phantom, coiled up, D — 1-dimensional 
polymeric network of size L± embedded in d — 1 dimensions, with the radius of 
gyration Rq{L±) ~ D ^ 2 . In the physical dimensions (D — 2 and d — 3) this in 

1/2 

particular gives a coiled up ideal polymer of length L± with Rq ~ L_/_ , as expected. 



Similar analysis for the tubule roughness h rms gives 
h 2 rms = 2(d-D) 



k B T f dq y I ^ 



k B T 



1 (1 _ Jly 1 

V 

ft „i 

y 



L?- 1 J L ? (27T) n{q y )q y 

(44) 



d D - 1 q±dq y (\-e tq y L 



In contrast to Rq, only the = zero mode (first term) and bulk modes contribute 
to h rmSl giving 

h rrns (L ±1 L y ) = L y S h {L y /Ll) (45) 
with, for phantom membranes, 

, 5 -2D 1 

^=—2-^-2- (46) 
and the asymptotics of Su{x) given by 

„ , , f constant for x — ► , . 

^i^/H, for x — * co (4?) 

Equations 1401 and 1451 give information about the tubule roughness for arbitrarily 
large size L±_ and L y , and arbitrary aspect ratio. For the physically relevant case of 
a square membrane L± ~ L y ~ L — > co, for which L y >> in contrast to Rq 
(where bulk modes dominates), qj^ = zero mode dominates, leading to 



L 



C+(D-l)/2z 

v_ 

L (D-l)/2 



OC L C+(D-l)(l-z)/2z ^ (4g) 



Equations IT51 HSl then give, for a D = 2 phantom tubule, £ = 1/2, z = 1/2 



L 



3/2 



V 



L 1 / 2 



(49) 



and therefore predicts for a square phantom membrane wild transverse tubule un- 
dulations 

Kms ~ L , (50) 
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consistent with simulation:^! that find h rrns ~ L 7 , with 7 = 0.895 ± 0.06. As with 
i?G, it would be interesting to test the full scaling law Eq03]by simulating non- 
square membranes, and testing for the independent scaling of h rms with L y and 
L±. Note that, unlike Rg, according to Eg 1491 h rms will show immediate growth 
(reduction) when one increases (decreases) L y at fixed L±. 

The above discussion also reveals that our earlier conclusions about the lower 
critical dimension D/ c for the existence of the tubule are strongly dependent on how 
L± and L y go to infinity relative to each other; i.e., on the membrane aspect ratio. 
The earlier conclusion that D; c — 3/2 only strictly applies when the bulk modes 
dominate the physics, which is the case for a very squat membrane, with L y w L z ± , 
in which case L y « L±. For the physically more relevant case of a square phantom 
membrane, from the discussion above, we find that tubule phase is just marginally 
stable with Di c = 2~ 

Equations 1401 and 1451 also correctly recover the limit of L z ± = constant <C L y — > 
00, where the tubule simply becomes a polymer (ribbon) of thickness Rg{L±) of 
length L y embedded in d — 1 dimensions. These equations then correctly recover 
the polymer limit giving 

h rms « Lp{L y /L P fl 2 , (51) 

with Lj_-dependent persistent length Lp(L±) oc _1 . So, as expected for a phan- 
tom tubule, if L± does not grow fast enough (e.g. remains constant), while L y — > 00, 
the tubule behaves as a linear polymer and crumples along its axis and the distinc- 
tion between the crumpled and tubule phases disappears. 



2.6. Self-avoidance in the tubule phase 

Self-avoidance is an important ingredient that must be included inside the tubule 
phase. Detailed analysis of self-avoidance overturns arguments in Sec l2.5.l1 and 
leads to anomalous elasticity in the bending rigidity modulus «(q). Self-avoidance 
therefore also modifies the values of other shape exponents, while leaving the scaling 
form of correlation functions unchanged. 

In the y-tubule phase the self- avoidance interaction F$a from EqJSJ reduces to 

F SA = v fdyd^x^-^S^ih^y)-^^)) , (52) 

with v — b/2£ y . 

2.6.1. Flory theory 

The effects of self-avoidance can be estimated by generalizing standard Flory 
argument^^l from polymer physiciPSI to the extended tubule geometry. The total 
self- avoidance energy scales as 

E SA « Vp 2 , (53) 
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where V ~ R d Q 1 L y is the volume in the embedding space occupied by the tubule 
and p = M/V is the embedding space density of the tubule. Using the fact that the 
tubule mass M scales like L^~ X L V , we see that 



L q ,L 



2(13-1) 



Esa « y > (54) 

Using the radius of gyration Rq oc L 1 ^, and considering, as required by the 
anisotropic scaling, a membrane with L± oc Ly, we find that Esa oc Ly SA around 
the phantom fixed point, with 

Asa = 1 + 4(D - 1) - (d - , (55) 

Self-avoidance is relevant when \$a > 0, which, from the above equation, hap- 
pens for v = i/ph = (5 — 2D) /A when the embedding dimension 

For D = 2-dimensional membranes, d^T = 11. Thus, self-avoidance is strongly 
relevant for the tubule phase in d = 3, in contrast to the flat phase. 

We can estimate the effect of the self- avoidance interactions on Rq in Flory 
theory, by balancing the estimate Eg 1541 for the self-avoidance energy with a similar 
estimate for the elastic energy: 

Eelastic = t (— — J 1 L y . (57) 

Equating E e i ast i C with Esa, we obtain a Flory estimate for the radius of gyration 
Rg- 

R g (L ± ) ck L v l , v F = , (58) 

which should be contrasted with the Flory estimate of Vp = (D + 2)/(d + 2) for the 
crumpled phase. For the physical case D = 2, Eg 1581 gives 

Rg oc L 3 / 4 , (59) 

a result that is known to be exact for the radius of gyration of a D — 1-polymer 
embedded in d — 2-dimensions@3 Since the cross-section of the D — 2-tubule, 
crudely speaking, traces out a crumpled polymer embedded in two dimensions (see 
FigOU, it is intriguing to conjecture that = 3/4 is also the exact result for the 
scaling of the thickness of the tubule. 



2.6.2. Renormalization group and scaling relations 

A new significant complexity that arises and is special to the tubule phase (as 
compared to the crumpled phase) is the simultaneous presence of local elastic 
and nonlocal (in the intrinsic space) self- avoidance nonlinearities. Above Flory 
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mean-field analysis (that ignores elastic nonlinearities) is nicely complemented by a 
renormalization group approach that can handle this complexity of the full model 
F = F e i + FgA, Eas l24l52l Although, as we argued above, elastic nonlinearities are 
irrelevant in a phantom tubule, in a physical self-avoiding tubule, they are indeed 
important and lead to an anomalous bending rigidity elasticity. 

The correct model, that incorporates the effects of both the self-avoiding inter- 
action and the anharmonic elasticity, is defined by the full tubule free energy 

2-i 



F = \i d D - l .r ±l hj 



K (d 2 h) 2 + t(d^h) 2 + g^uf + g y (dyu + \{d y hf 

+ v fdyd D - 1 X ± d D - 1 x' ± 5^- 1 \h{x±,y)-h(x' ± ,y)) . (60) 

To assess the role of elastic (g v ) and self-avoiding (v) nonlinearities we implement 

48 

standard momentum-shell RGI-^We integrate out perturbatively in g y and v short- 
scale fluctuations of modes u(q) and h(q) within a cylindrical shell Ae~ ; < q± < A, 
— oo < q y < oo, and anisotropically rescale lengths (xj_, y) and fields (/i(x), u(x)), 
so as to restore the ultraviolet cutoff to A: 

x±=e'x' ± , y = e zL y', (61) 
h{x) = e' A h\x!) , u(x) = e (2,y - 2 V(x') , (62) 

where we have chosen the convenient (but not necessary) rescaling of the phonon 
field u so as to preserve the form of the rotation-invariant strain operator d y u + 
^(d y h) 2 . Under this transformation the free energy returns back to its form, Ea l60l 
but with effective length-scale (I = logLx) dependent coupling constants deter- 
mined by 



dt 
dl 
dn 
~dl 
dg y 
dl 
dg± 
dl 
dv 

~dl 



[2v + z + D-3-f t (v)]t, (63) 

[2v-3z + D-l + f K {g y ,g x )] K , (64) 

[4v-3z + D-l-f g (g y )]g y , (65) 

[Av-z + D- 3}g x , (66) 

[2D - 2 + z - (d- \)v- f v (v)}v , (67) 



where the various /-functions represent the graphical (i.e., perturbative) corrections. 
Since the self-avoiding interaction only involves h, and (for convenience) the param- 
eters in the h propagator (t and k) are going to be held fixed at 1, the graphical 
corrections coming from self-avoiding interaction alone depend only on the strength 
v of the self-avoiding interaction. Therefore, to all orders in v, and leading order 
in g y , f t (v) and f v (v) are only functions of v and f K (g v ,gx) and f g {g y ) are only 
functions of g y and g±_ . 
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It is important to note that g± suffers no graphical corrections, i.e., Eg 1661 is 
exact. This is enforced by an exact symmetry 

u(x±,y) u(x.±,y) +x( x -l) ) (68) 

where x( x -0 is an arbitrary function of xj_, under which the nonlincaritics in F are 
invariant. 

We further note that there is an additional tubule "gauge" -like symmetry for 
9y = 

h(x x ,y)^h(xL X ,y) + $(y) , (69) 

under which the only remaining nonlinearity, the self-avoiding interaction, being 
local in y, is invariant. This "tubule gauge" symmetry demands that f K (9y = 
0,g±) = 0, which implies that if g y = 0, there is no divergent renormalization 
of k, exactly, i.e., the self-avoiding interaction alone cannot renormalize k. This 
noTi-renormalization of k by the self-avoiding interaction, in a truncated (unphysi- 
cal) membrane model with g y = , ha s been verified to all orders in a perturbative 
renormalization group calculation!^ 

To see that v and z obtained as fixed point solutions of Eos lfi 31671 have the same 
physical significance as the v and z defined in the scaling expressions EqsHfor the 
radius of gyration R G and tubule wigglyness h rms , we use RG transformation to 
relate these quantities in the unrenormalized system to those in the renormalized 
one. This gives, for instance, for the radius of gyration 

R G {L ± ,L y ; t,K,...) = e" l R G {e~ l L ± ,e- zl L y ; t(l), K (l), . . .) , 

(70) 

Choosing I = /* = logLj_ this becomes: 

Rg(L±, L y ; *,«,...) = L v ± R G (l,L y /L z ± ; t{h), K (h), . . .) . (71) 

This relation holds for any choice of the (after all, arbitrary) rescaling exponents 
v and z. However, if we make the special choice such that Eas l63l67l lead to fixed 
points, ... in Ea l71l go to constants, independent of I* (and hence L±), 

as L± and hence I*, go to infinity. Thus, in this limit, we obtain from Eg 1711 

R G (L ± ,L y ; t,K,...)= LlR G (l, L y /L z ± ;U, «.,...), (72) 

where are the fixed point values of coupling constants. This result clearly 

agrees with the scaling forms for R G , Eq^ (with analogous derivation for h rms ) 
with the scaling function given by Sr{x) = R G (1, x; t*, n*,g*,v*). 

The recursion relations Eas l63l(37l rcproduce all of our earlier phantom membrane 
results (when v = 0, leading to /„ = 0), as well as the upper-critical embedding 
dimension = (6D — l)/(5 — 2D) for self-avoidance predicted by Flory theory, 
Eg 1551 and the upper critical intrinsic dimension D uc = 5/2 for anomalous elastic- 
ity for phantom membranes (determined by eigenvalues of dimensionless couplings 
corresponding to v and g y ). They also reproduce all of the Flory theory exponents 
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under the approximation that graphical corrections to t and v are the same, i.e., 
that ft(v*) = fv(v*). 

To analyze the renormalization of k in a self-avoiding tubule, we focus once 
again on the nonlocal interaction the non-local interaction F an h, Ea l27l mediated 
by integrated out phonons, with a kernel 

^(<0 = i ¥^4 2 > (73) 



9±Q 



± 



whose long-scale scaling determines renormalization of n. If g y (<l)qy >> .9±(q)<zi 
(as we saw for a phantom tubule) then at long scales 14 (q) « g±q±/qy ~ in 
the relevant limit of q± ~ q^. Simple power counting around the Gaussian fixed 
point then shows that this elastic nonlinearity is irrelevant for D > D uc = 3/2, and 
therefore /* = r\ K — for a physical D = 2-dimensional tubule, as we argued in 
Secl231 

On the other hand, if the scaling is such that g±(q)q? L dominates over ff v (q)g?, 
then 14 (q) ~ 9y, i-e. a constant at long length scales. Simple power-counting then 
shows that this coupling is relevant for D < D uc = 5/2 and the bending rigidity 
modulus of a D = 2-dimensional tubule is anomalous in this case. 

In the RG language, the relevance of 14 is decided by the sign of the renormal- 
ization group flow eigenvalue of g±(l) in Eg 1661 

X g± = 4u - z + D - 3 , (74) 

which is exactly determined by the fixed-point values of v and z, since g± suffers 
no graphical renormalization. 

As discussed in previous sections, for a phantom tubule (or for d > d uc ) v = 
(5 - 2D) /A and 2 = 1/2. For d below but close to d^ = (6D - l)/(5 - 2D) (= 11 
for D = 2), these values are modified by the self-avoiding interaction, but only by 
order e = d — d^ . Hence a D — 2-dimensional tubule, embedded in d dimensions 
close to <f^g =llj A gx w — 1/2 and g±(l) flows according to 

^ = l-\ + 0(e)] 9± , (75) 

i.e. g± is irrelevant near d — 11 (for e -C 1), 14(q) ~ 9±q±/qy ~ <?y °^' ) i s irrelevant 
to the bend elasticity for a D = 2-dimcnsional tubule embedded in high dimensions 
near df^ = 11, and, hence, f K — r\ K in Ea 1641 vanishes as the fixed point. Therefore, 
in these high embedding dimensions the full model of a self-avoiding tubule reduces 
to the linear elastic truncated model with self-avoiding interaction as the only non- 
linearity, that can be nicely studied by standard expansion in e = d^ — As 
we discussed above, the "tubule gauge" symmetry guarantees that in this case the 
self-avoiding interaction alone cannot renormalize k, i.e., f K = r\ K = for d near 
d^. This together with Eq[M]at its fixed point leads to an exact exponent relation 



z=^(2v + D-\) , 



(76) 
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valid for a finite range < d < df^ of embedding dimensions, and for phantom 
tubules in any embedding dimension. Where valid, it therefore reduces the tubule 
problem to a single independent shape exponent. 

However, this simple scenario, and, in particular, the scaling relation Eg 1761 is 
guaranteed to break down as d is reduced. The reason for this is that, as d decreases, 
v increases, and eventually becomes so large that the eigenvalue X g± of g± changes 
sign and becomes positive. As discussed earlier, once this happens, the nonlinear 
vertex Eg 1731 becomes relevant, and k acquires a divergent renormalization, i.e., 
fn 7^ 0; an d bend tubule elasticity becomes anomalous. 

Now, it is easy to show, using Eo l76l and a rigorous lower bound v > (D — 
l)/(d — 1) inside X g± , that it must become positive for d > d l ^{D) with 

d * {D) = T^F = 7/2 ' f° r D = 2. (77) 

Hence, for the case of interest D = 2 critical dimension is bounded by 7/2 from 
below. In fact, Flory and e = 11 — d estimates indicate that d*(D = 2) w 6p21 

We therefore conclude that in a physical D — 2-dimensional self-avoiding tubule, 
embedded in d = 3 < d* « 6, anharmonic elasticity F an h is important at long scales 
and leads to anomalous and divergent bending rigidity /t(q) and <?j,(q) 

g y (q) = qfS g {q v /ql) , «(q) - q-"S K (q v /q> ± ) , (78) 

with 

zVk = fn{g* y ,9*x) , zr] u = f g (g*) . (79) 

and asymptotic forms of scaling functions 

S g (x -> 0) -» as-"" , S K (x -> 0) -» . (80) 

Another consequence is the breakdown of the exponent relation, Eg 1761 that is 
replaced by two exact relations holding between four independent exponents z, v, 
rj K , and rj u 

z = -^—{2v + D-l), z =-±—(4 V + D-l), (81) 

which automatically contain the rotational symmetry Ward identity 

2r? K + Vu = 3 - (D - l)/z , (82) 

formally arising from the requirement that graphical corrections do not change the 
form of the rotationally invariant strain operator d y u + ^(d y h) 2 . The existence of a 
nontrivial > 3 and its consequences are summarized by Figs l6l7l 

The physics behind above somewhat formal derivation of exponent relations 
(Ward identities) can be further exposed through a simple physical shell argument. 
As can be seen from FiglHI bending of a tubule of radius Rg into an arc of radius 
R c induces an in-plane strain energy density g y (L y , L±)(Rg(L v ) / R c ) 2 ■ Interpreting 
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Fig. 6. Schematic illustration of change in relevance of g±(l) at d*. For embedding dimensions 
below d* (which includes the physical case of d = 3), <?j_(0 becomes relevant and (among other 



phenomena) leads to anomalous bending elasticity with re(q) 1 
scales. 



ly 



that diverges at long length 



this additional energy as an effective bending energy density n y (L±, L y )/R^, leads 
to the effective bending modulus K y (L±, L y ), 

Ky(L±,L y ) ~ g y (Lj_, L y )Rc(L±, L y ) 2 . (83) 

This leads to a relation between the scaling exponents 

2u = z{r] K + r lu ) . (84) 

that is contained in Ecisl81l obtained through RG analysis above. 

We note, finally, that all of the exponents must show a jump discontinuity at <i*. 
Therefore, unfortunately, a n ex trapolation from e = 11 — d-expansion in a truncated 
model with linear elasticity^ down to the physical dimension of d — 3 (which is 
rigorously below d*) gives little information about the properties of a real tubule. 

The computations for a physical tubule must be performed for d < d*, where 
both the self-avoidance and the elastic nonlincaritics are both relevant and must be 
handled simultaneously. As we discussed above, for d < d*, the eigenvalue X g± > 0, 
leading to the flow of g±(l) to infinity, which in turn leads to V/,.(q) — g y . Physically 
this regime of g± — > oo corresponds to freezing out the phonons u, i.e. setting u = 
in the free energy F[h,u], EqEOI with the resulting effective free energy functional 
for a physical self-avoiding tubule given by 



x±dy 



K{dlhf+t{d^hf 



dyd D 1 x^d 



>-'- 1 x' ± S( d - 1 \h(x ± ,y)-h( X ' ± ,y)) 



(85) 
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Fig. 7. Schematic of the tubule "phase" diagram in the embedding d vs intrinsic D dimensions. 
Self-avoiding interaction becomes relevant for d < d^(D) = (GD — 1)/(5— 2D), (=11, for D = 2). 
Below the d t (D) curve (for which the lower bound is d l ^(D) = (AD — l)/(4 — D)) the anharmonic 
elasticity becomes relevant, leading to anomalous elasticity with a divergent bending rigidity. 




Fig. 8. Illustration of the physical mechanism for the enhancement of the bending rigidity k by 
the shear g y elasticity. To bend a polymerized tubule of thickness Rq into an arc of radius R c 
requires Rq / R c fraction of bond stretching and therefore costs elastic shear energy, which when 
interpreted as bending energy leads to a length-scale dependent renormalization of the bending 
rigidity k and to the Ward identity Eg 1841 
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Since one must perturb in g y around a nontrivial, strong coupling fixed point with 
v* = O(l), unfortunately, no controlled analysis of above model exists todate and 
remains a challenging open problem. Nevertheless, above RG analysis combined 
with Flory estimates and exact exponent relations provides considerable information 
about shape fluctuations and elasticity of a polymerized tubule. 

2.7. Phase transitions 

Now that we have solidly established the properties of the three phases of anisotropic 
polymerized membranes, we turn to analysis of phase transitions between them. As 
discussed in the Introduction, a direct crumpled-to-flat transition is highly non- 
generic for anisotropic membranes, as it has to be finely tuned to the tetracritical 
point illustrated in Fig|2 If so tuned this transition will be in the same universality 
class as for isotropic membranes (where it is generic) 

Here we will focus on the new crumpled-to-tubule and tubule-to-flat transi- 
tions, which are generic for membranes with any amount of anisotropy. As for the 
tubule phase itself and the crumpled-to-flat transition, a complete analysis of these 
transitions must include both elastic and self-avoiding nonlincarities, a highly non- 
trivial open problem. Below we will instead present an incomplete solution. First 
we will present an RG analysis of a phantom (non-selfavoiding) membrane, focusing 
on the much simpler crumpled-to-tubule transition. We will then complement this 
study with a scaling theory and Flory approximation of the crumpled-to-tubule and 
tubule-to-flat transitions in a more realistic self-avoiding membrane. 

2.7.1. Renormalization group analysis of crumpled-to-tubule transition 

The crumpled-to-(y-)tubule (CT) transition takes place when t y — > 0, while t± 
remains positive. This implies that the CT critical point is characterized by highly 
anisotropic scaling q± cx q y . It leads to a considerable simplification of the full free 
energy defined in Eq|S]down to 

= \J d D - l x ± dy n y (d 2 y r) 2 + t x {d^ff + t y (8 y r) 2 + (8 y r ■ d y rf 

(86) 

The standard 0{d) 4> 4 mo del form of F facilitates systematic analysis using con- 
ventional RG methodsPSI Because of the strong scaling anisotropy of the quadratic 
pieces of the free energy, it is convenient to rescale lengths xj_ and y anisotropically: 

x±=x' ± e', y = y'e zl 1 (87) 

and rescale the "fields" r (x) according to 

r(x) = e*¥(x') • (88) 

Under this transformation 

Ky(l) = Kye^- 1 - 3 ^ 1 , tx(l) = txe^ D - 3+z+2 ^ 1 . (89) 
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Requiring that both K y and t± remain fixed under this rescaling (zeroth order RG 
transformation) fixes the "anisotropy" exponent z and the "roughness" exponent 
X- 

z = \, x=(5/2-£>)/2. (90) 

Although this choice keeps the quadratic in r part of F, Ea l86l unchanged, it does 
change the quartic piece: 

u yy(l) = u yy e ^ + * = u yy e ^ ^ ^ ! (91) 

and shows that below the upper critical dimension D uc = 5/2, the Gaussian critical 
point is unstable to elastic nonlinearities, that become comparable to harmonic 
elastic energies on scales longer than the characteristic length scale Uf_ 

Lf=(^L) . (92) 

To describe the new behavior that prevails on even longer length scales requires a 
full-blown RG treatment. 

Such an analysis^' leads to corrections to the simple rescaling of K y , t±, and t y 
(due to u yy non-linearity), characterized by "anomalous" exponents rj K , rj t , and 69: 

Ky(l) = Kye^- 1 - 3 ^^ 1 , (93) 
t ± (l) = t±e (D-3+z+ Vt +2 X )i 5 ( 94) 

ty (l) = t y e( D - l - z - se+2 x» = t y e Xtl , (95) 



that give 

_ 2 - rjt = 10 - 4D + r] K {D - 3 + f] t ) -3^ 
4 - Vk 1 X 8-2 Vr 

valid at the new nontrivial CT critical point. 

Once the values of rjt, r] K and x a t the critical point are determined, the renor- 
malization group gives a relation between correlation functions at or near criticality 
(small t y ) and at small wavectors (functions that are difficult to compute, because 
direct perturbation theory is divergent) to the same correlation functions away from 
criticality and at large wavevectors (functions that can be accurately computed us- 
ing perturbation theory). For example the behavior of the correlation lengths near 
the transition can be deduced in this way: 

U(t v ) = e l U(t y e Xtl ) , ty(t y ) = e*%{t y e x * 1 ) , (96) 

assuming that a critical fixed point exists and all other coupling constants have 
well-defined values at this point. Taking t y e Xtl w 1 then gives 

U(ty)^at~^, iy{t v )nat- v y , (97) 
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where a ~ £(1) is the microscopic cutoff and, 

" ± = h = 2(2-^-^-^(2-^-^)' ^ = ^ (98) 

are correlation length exponents perpendicular and along the tubule axis. 

The anomalous exponents can be computed by integrating out short-scale de- 
grees of freedom perturbatively in u yy . Standard analysis shows that indeed there is 
a nontrivial critical point (at a finite value of u* y = 0(e), e = 5/2 — D)), at which, 
to all orders rj t = 0, and to one-loop order, for a physical membrane (£) = 2, d = 3, 
i.e., e= 1/2) 

= 0, z = 1/2, x = 1/4, i>l « 1-227 , v y « 0.614 . (99) 

It is interesting to note that, in contrast to the treatment of crumpled-to-flat tran- 
sition in isotropic membraneJ^El where the critical point was only stable for an 
unphysically large value of the embedding dimension d > 219, the critical point 
characterizing the crumpled-to-tubule transition discussed here is stable for all d. 
Furthermore, the relatively small value of e = 1/2 (in contrast to for example e = 2 
for flat phase and crumpling transition) gives some confidence that above critical ex- 
ponents for the CT transition in a phantom membrane might even be quantitatively 
trustworthy. 



2.7.2. Scaling theory of crumpled-to-tubule and tubule-to-flat transitions 

Above RG analysis for the CT transition in phantom membranes, can be nicely com- 
plemented by a scaling theory combined with Flory estimates, that can incorporate 
both the elastic and self-avoiding nonlinearities, as we now describe. 

Standard scaling arguments, supported by RG analysis and matching calcula- 
tions (see e.g., Sec l2.5|) suggest that near the crumpled-to-tubule transition, for a 
square membrane of internal size L, membrane extensions R y and Rg along and 
orthogonal to the developing tubule axis should exhibit scaling form: 

Rg, v = L^ v f G ,v{tyL+) , 

y > 0,L » £ ct 

L«U (100) 
t y <0,L » £ ct 

where subscripts t, c and ct refer to tubule, crumpled and tubule-to-crumpled tran- 
sition, respectively, and £ c j cx |i y | -1 ^ is a correlation length for the crumpled-to- 
tubule transition, t y = (T — T ct )/T ct , T ct is the crumpled-to-tubule transition tem- 
perature, with t y > corresponding to the crumpled phase. Consistency demands 
that exponents 7^y^_ are given by 

,.G,V ,G,y G,y 

7 G,y = Vc-V ct = V t -Vet ^ 
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The asymptotic forms in Ea llUUI are dictated by general defining properties of the 
phases and the CT transition. For example, scaling of both R y and Rq like L v " , with 
the same exponent v c in the crumpled phase is rooted in the isotropy of that phase. 
In contrast, extended and highly anisotropic nature of the tubule phase dictates 
that f t G ^ v\ = 1. The anisotropy is, however, manifested even in the crumpled 
phase through the different temperature-dependent amplitudes of Rq and R y , with 
the aspect ratio R y /Rc actually diverging as T — ► T^, and membrane begins to 
extend into a tubule configuration. The former of these vanishes as t y — ► + (since 
the radius of gyration in the tubule phase is much less than that in the crumpled 
phase, since v t < v c ), which implies 7^ > 0, while the latter diverges as t y — ► + , 
since the tubule ultimately extends in that direction, which implies 7+ < 0. 

We will now show how these general expectations are born out by the Flory 
theory. Following analysis very similar to that done in Sec l2.6.il but keeping track 
of temperature-dependent order parameter Q y , we find that the Flory approximation 
to the free energy density in Eg 1861 supplemented with self-avoidance is given by 



fFl[RG,Cy]=tyCy+u vy Ct+tx[Y^) +v^ 1 ^rT- (102) 



Minimizing this over Rg, gives 



\t±C 



i/(d+i) 



with the tubule exponent v t — found earlier, but now with additional tem- 
perature and L± dependence of Rq through ( y (t y ,L±) that interpolates between 
tubule, crumpled and critical behavior. Inserting this expression for Rq into Ea ll()2l 
gives 



d+T 2(d-D) 

f-2 



fFl [ty] = tyCy + U yy Cy + tf 1 [± j L ± d+1 . ( 1 0-1 ) 

Minimizing fpi [£ y ] in the crumpled phase (t y > 0) gives 



2 d _l \ 2P+2J 



Cy « ) L ± d+2 , (105) 

that, as expected vanishes in the thermodynamic limit. For a square (L x L) mem- 
brane this then gives for R y = C, y L y and Rq (using Ea ll03|) 

RyOcty^Lj^ , Raocty^Lj^ , (106) 

which, after comparing with the general form, Eg 1 1001 gives 

D+2 y _ d+1 G _ 1 
Vc ~ ~d + 2 ' 7+ ~ ~2(d + 2) ' 7 + ~ 2(d + 2) ' ( ] 
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v c reassuringly agrees with the well-known Flory result for the radius of gyration 
exponent v c for a D-dimensional manifold, embedded in d dimensions EE2E3I23I 
and 7^.' G special to crumpled anisotropic membranes. 

As anticipated earlier, 7^ ^ 7+ implies that intrinsically anisotropic membrane 
are qualitatively distinct from isotropic ones even in the crumpled phase, as they 
exhibit a ratio of major to minor moment of inertia eigenvalues (related to Rq/R v ) 
that diverges as CT transition is approached. 

Now for the tubule phase, characterized by t y < and a finite order parameter 
C, y > 0, last term in fpi, Ea ll(J4l is clearly negligible for L± > £ cr and simple mini- 
mization gives C, v ~ \/\ty\/u y y, which, when then inserting into R v ,g and comparing 
with the general scaling forms gives for a square membrane 

^-'•^--l^-TTT'^—wnr)- (108 > 

Finally, right at the crumpled-to-tubule transition, t y — 0, minimization of 
fFi[( y ] gives 

t v « L"^ (109) 

which, when inserted in R Vt G gives the advertised critical scaling forms with expo- 
nents 

„ D + d + 3 G 2D + 3 2(d-D) 

^ = ^2dT^> ' * = ^dT3T- (no) 

that are reassuringly consistent with our independent calculations of exponents 
7_l_'* j Vc-, v t ' V 1 anf i v ct' V using exact exponent relations above. For the physical case 
of a two dimensional membrane embedded in a three dimensions, (D = 2, d = 3) 

^ c = 4/5, ^ = 7/9, 1^ = 8/9, i* = 3/4, (111) 
7 G = 1/10, 7* = -2/5, 7 G = -l/8, 7- = V2, = 2/9, 

The singular parts of other thermodynamic variables obey scaling laws similar 
to that for Ra,y, Ea llUUI For example, the singular part of the specific heat per 
particle C v ~ 777 Jjt (h^y^y^ D 2 )^ usm S Ea llOOl exhibits scaling form 



a 




ty > 0,L >> tct 

L«U (112) 

ty <Q,L» £,Ct 



with ^^[/'(x)], and 



P = 2i& - 2 + cj> , a± = -2j y ± + 1 , (113) 

rv = 

d + 2 



13 = , a+= M + 3 , a_ = , Flory theory. (114) 
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This leads to the unusual feature that outside the critical regime (i.e. for L >> 
£ c t), the singular part of the specific heat above the crumpled-to-tubule transition 
vanishes in the thermodynamic limit like L~ Q +^ ~ L _2 ( d_ - D )/( d + 2 ) ~ L~ 2 / 5 . Similar 
behavior was also found for the direct crumpled-to-flat transition by Paczuski et 
al.P 

We now turn to the tubule-to-flat (TF) transition. On both sides of this transi- 
tion, R y = L y x 0(1). Therefore only Rq exhibits critical behavior, which can be 
summarized by the scaling law: 

R G = L^f V f ± (t x L^f) , 

t ± >0,L»£ tf 
oc < L Vt f , L « £ tf (115) 

[\t±\"'-L,t±<0,L»{t / 

where t± = (T — T t f)/T t f, tj_ > corresponds to the tubule phase, £tf oc \t±_\~ 1 ^ tI 
is the correlation length for this transition, and the exponents obey the scaling 
relations 

ii = iyt - f*/)M/, il = (i - «v)M/ ■ (H6) 

In Flory theory we find: 

&/ = l/3, v tf = 5/6, 7 ^ = -1/4, 7 l / = -l/2. (H7) 
The singular part of the specific heat again obeys a scaling law: 

C v = L 2 " t f +<l > t f- 2 g(t ± L^f) , 

ftI a +L-"+, t x >0,L»€ t/ 

OC ^ L 2, t/ +0 t/ -2 ; L <<; £ {/ (11 8 ) 

where, in Flory theory, 

a^=3/2, a£f=0, k + = 1/2, 2i/ t / + 2 = k_ = 0. (119) 

Thus, again, the singular part of the specific heat vanishes (now like £ -1 ' 2 ) in 
the thermodynamic limit above (i.e., on the tubule side of) the transition, while it 
is 0(1) and smooth as a function of temperature in both the critical regime and in 
the flat phase. 

3. Random Heterogeneity in Polymerized Membranes 
3.1. Motivation 

Soon after a general picture of idealized homogeneous membranes was established, 
much of the attention turned to effects of random heterogeneity on conformational 
properties of polymerized membranes !^ As with many other condensed matter 
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systems (e.g., random magnets, pinned charged density waves, vortex lattices in 
superconductors a main motivation is that random inhomogeneity is an in- 
evitable feature in most physical membrane realizations. As illustrated a cartoon 
of cellular membrane wall (Fig. 1.4 in lectures by Stan Leibler), functional proteins, 
nanopores (controlling ionic flow through membrane), and other heterogeneities 
(with sincere apologies to biologists for such crude physicist's terminology) are in- 
corporated into a cellular lipid bilayer. Holes or tears, random variation in the local 
coordination number (disclinations), dislocations, grain-boundaries, and impurities 
incorporated into fishnet-like biopolymer spectrin network attached to cellular wall 
(as e.g., in red blood cells) are also substantial sources of inhomogeneity in mem- 
branes. Such defects in the two-dimensional polymer network are also inadvertently 
introduced during photo-polymerization of synthetic polymerized membranes. If one 
is only interested in statistical conformational properties of such membranes, these 
frozen-in heterogeneities can be treated as random quenc hed disorder, similar to 
treatment of impurities in other condense d m atter contexts!^ 

Consistent with theoretical predictions the fact that quenched internal dis- 
order can drastically modify the conformational thermodynamics of polymer- 
ized membranes was first demonstrated experimentally by Mutz, Bensimon and 
Brienne^l They observed that partially (heterogeneous) polymerized vesicles un- 
dergo upon cooling a "wrinkling" transition to a folded rigid glassy stucture that re- 
sembles a raisin. Natural interpretation of this important experiment as an evidence 
for a transition towards a crumpled spin-glass-like state provided strong motivation 
for further theoretical studies of quenched disorder in polymerized membranes. 

Below we will describe a generalized model that includes effects of quenched 
disorder in a polymerized membrane and will show that (as in-plane anisotropy, 
discussed in previous section) it has drastic qualitative effects on long-length con- 
formational properties of a polymerized membrane. 



3.2. Model of a heterogeneous polymerized membrane 

It is clear that above sources of heterogeneity lead to local random in-plane di- 
lations and compressions and can therefore be modelled by local random stresses 
— a a fl(ii)u a p. Geometrically this can be understood as a random preferred back- 
ground metric g°p(x) = 5 a/3 + ^(x) with strain u a p = \(g a f3 - 3a/3( x )) "mea- 
sured relative to this deformed state, and metric g a p = d a f ■ dpr seeking to relax 
to <7°^(x)). These local in-plane stresses can and will be relaxed by buckling of 
the membrane into the third dimension that will tend to screen the elastic inter- 
action, thereby lowering the elastic energy by partial trade-off between in-plane 
clastic energy and membrane bending (curvature k) energy. However, as illustrated 
in Fig|5Ja), because such randomness respects the reflection symmetry relating two 
sides of the membrane, the induced puckering will locally break this Ising symmetry 
spontaneously, and in a way specific to each configuration of disorder. 

A qualitatively distinct form of quenched disorder that explicitly breaks reflection 
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(a) 




Fig. 9. A cartoon of a bilayer membrane with a reflection- (a) symmetric and (b) asymmetric 
inclusions that can be modelled by two qualitatively distinct types of disorder, the random stress 
and random mean curvature disorder, respectively 



symmetry arises from asymmetric inclusions of the type illustrated in FigEfb). 
These lead to a local preferred extrinsic curvature, that, in the flat phase is described 
by ~V 2 h ■ c(x). 

Membrane defects will also of course lead to heterogeneous elastic moduli, k(x), 
/x(x), and A(x). However, it can be shown that such weak heterogeneity (i.e., as long 
as the average value of these elastic constants remains larger than their variance, 
that is they are predominantly positive) has no qualitative effects on membrane 
long-scale conformations. Consequently, all effects of membrane heterogeneity can 
be modelled by just two types of quenched disorder, random stress and random 
curvature. This is perhaps not surprising given the aforementioned analogy of mem- 
branes with ferromagnet (with membrane normal n playing the role of a spin S), 
where too, random bond (that respects S — > —S) and random field (that is odd un- 
der S — > — S) are the only two qualitatively important types of quenched disorder. 
One qualitatively important distinction from random magnets that we can already 
anticipate at this point is that the curvature disorder, that couples to the gradient of 
the order parameter h is far weaker perturbation than its ferromagnetic analog, the 
random field disorder that couples directly to magnetization. This distinction will 
lead to importance of the curvature disorder below D uc = 4, (same as the random 
stress disorder and therefore competing with it) contrasting with the upper-critical 
dimension of D uc — 6 of the random-field in a ferromagnetic 

The general model of heterogeneous membrane is therefore described by an 
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effective Hamiltonian F[u a , h] 



F[u a ,h] = J d D :< 



-(V 2 /i — ) + m(«q/3) 2 + 2 {u aa ) 2 - fJ>u a 0T) a p(x) - -up/3r) aa (x) 



c(x) \ 2 . * 9 A . , 9 , . A 

. . / \ .- . i .~ .... ...... _ , ( 

(120) 

where quenched disorder fields c(x) and r] a p{yi) can be characterized by zero-mean, 
Gaussian statistics with second moment given by: 



Ci(x) Cj (0) = «SyA c (x) , (121) 

1 



77 Q/3 (x)?7 75 (0) = (Ai(x) - — A 2 (x))S a pS 7 s + ^A 2 (yi)(S a7 Sp S + S ag S^) . (122) 

Another useful form of this model is obtained after the phonons u a are integrated 
out, which, at long length scales, can be done exactly since the phonons (unlike 
h) are not soft and therefore can be approximated by a harmonic elasticity. The 
resulting Hamiltonian is given by 

F[h] = J d D x [|(V 2 K - ^) 2 + \d a h ■ d h {2*1*1% + J^i$J$} d y h ■ dsh 



-iry^(x) (2^Pj 5 + J^-P^pr | dj h ■ dsh 



(123) 



where Pj^ = 8 a p ^r- and P^ = are transverse and longitudinal projection 

operators. For D = 2, d = 3 membrane F[h] simplifies considerably to: 



F[h] = J d 2 x 



(124) 



For a generic configuration of impurity disorder, the ground state is highly nontrivial 
as it is determined by simultaneous, but generically conflicting, minimization of the 
extrinsic and Gaussian curvature (R = i((V 2 /i) 2 — (d a dph) 2 )) terms 

V 2 h = i c (x) , P^d a hd p h = Pj^^x) . (125) 

Hi 

Long-scale properties of such ground state are amenable to statistical treatment, 
utilizing standard field theoretic machinery. 

3.3. Weak quenched disorder: "flat-glass" 
3.3.1. Short-range disorder 

For many (but not all; see below) realizations of heterogeneity discussed above, 
such as for example random membrane inclusions, the disorder is short-ranged, 
and therefore can be characterized by ^-function correlated disorder with variances 
A c (x) = A c <5 D (x), Ai, 2 (x) = A 1>2 ^(x). 

To understand the effects of quenched disorder it is helpful to first study the 
stability of the flat phase (described by the anomalous elastic fixed point, studied 
in Sec l2.4l) by performing a simple perturbative calculation in disorder and elastic 
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nonlinearities directly for a physical membrane (D = 2, d = 3). Standard analysis^! 
then leads to disorder and thermally renormalized bending rigidity k^: 

d?p K[q a P^(p)q ] 2 



(2tt) 2 K 2 |q + p| 4 
(2tt) 2 4 K |q + p| 



•(Ai+A,)/^ - rr.T- r . ( i26 ) 



The first, temperature-dependent correction that enhances k is identical to that 
of a homogeneous membrane and is responsible for the stability of the flat phase 
of polymerized disorder-free membranes.^ At low temperature the temperature- 
independent contributions dominate. The las t, r andom stress contribution leads to 
a divergent softening of the bending rigidit y!^ while the random curvature term 
works to increase the bending rigidit y I I, and thereby works to stabilize the flat 
phase through the "order-from-disorder" mechanism that is the zero-temperature 
analog of the thermal one discussed in Sec l2.4l 

Weak disorder should not affect the asymptotic behavior of membranes in the 
flat phase at sufficiently high temperatures, despite its importance at T = 0. To see 
this, assume the disorder is so weak that we can replace the elastic constants on 
the right hand side of Ea lT2*?jl bv wave- vector-dependent quantities kr(p) and Kn(p) 
renormalized only by thermal fluctuations in the way controlled by the disorder-free 
flat phase fixed point. As discussed in Sec l2.4l the se are expected to be singular at 
long scales, with kb.(p) ~ p^ v " and Kr{p) ~ p Vu P I 'H^ EI The expression for «Jj(?), 
Eci ll26l becomes 



d^p K R {p)[q a Pl {p)qp] 



> 2 4(q + p)|q + p| 4 



- {Al +A2) JW? 4« fl (q|q + p|4 ' (12?) 
= K R (q) [1 + const.A c g''" - const. (Ax + A 2 )q^} , (128) 

where we made use of the exact 2D exponent relatior P l ^ l ^ l 2rjf + -q u = 2 (a conse- 
quence of rotational invariance). Since rj K and r\ u are positive, at finite temperature 
weak quenched disorder just gives a subdominant nonanalytic correction to disorder- 
free result for Kn(q). Physically this finite temperature irrelevance of disorder comes 
from singularly soft in-plane elastic moduli and divergent bending rigidity, that, 
respectively facilitate screening of the in-plane stress induced by impurities and 
suppress wrinkling effects of curvature disorder. 

These perturbative arguments are supported by detailed renormalization group 
calculations controlled by an e = 4 — D-expansion that show that at finite temper- 
ature the Aronovitz-Lubensky fixed point is stable to weak quenched disorder!^! 
This is summarized by the RG flow equation of the coupling constants illustrated 
in FigUni 
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Fig. 10. RG flow diagram showing the irrelevance of random strain disorder (disorder variance 
scales to zero at long scales) near the disorder-free finite temperature fixed point P4 control- 
ling properties of the flat phase, and disorder induced instability of the flat phase at vanishing 
temperature (/^ e // °c pT). 



However, as is clear from the flow diagram, for sufficiently low-temperatures, 
even weak bare strain disorder becomes strong (A's flow to large values), invalidat- 
ing above perturbative argument. In this case a full low temperature RG analysis 
is necessary. As first discovered by Morse, Lubensky and Grest, it show d™^ 1 ^'" 
that interplay of random stress and curvature disorder leads to a new stable zero- 
temperature fixed point that controls long-scale properties of the disorder-roughened 
polymerized membrane. Similar to the thermally rough flat phase described by the 
AL fixed point, the resulting T — phase is characterized by a power-law roughness 
(with C < 1) about on average flat configuration. It therefore has all the ingredients 
of the "flat-glass" phase anticipated by Nelson and RadzihovskyP^I 

These RG results can be nicely complemented by the SCSA analysis.^ This 
can be done most effectively by first applying the replica formalism^!, that allows 
one to work with a translationally invariant effective Hamiltonian. To do this, one 
introduces n copies of fields h a and labeled by the replica index a, and integrates 
out the quenched disorder, thereby obtaining a replicated Hamiltonian. Assuming 
commutability of the thermodynamic and the n — > limits, the relation to the 
disorder-averaged free energy is established through the identity 



In Z = lim 



(129) 



where Z the partition function. 

The membrane roughness is characterized by the full disorder-averaged height 
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correlation function (that can also be related to replicated ones): 

<(Mx) - Mo)) 2 ) = ((h(x)-K(o))2) conn + (K(x)-Mo)) 2 , (130) 

~ A c \^ + A\^' , (131) 



where, respectively, the first (connected) and second contributions characterize 
thermal- and disorder-generated roughness, with corresponding roughness expo- 
nents £ and and the overbar denotes configurational (disorder) averages. The 
related exponents characterizing the Fourier transform of these parts of the height 
correlation functions are given by £ = (4 — D — rj R )/2 and C' = (4 — D — T)' K )/2. 
Analysis very similar to that done for homogeneous membranes in Sec l2.4.2l but 
with additional replica matrix structure leads t o a zero-temperature fixed point, 
that is marginally unstable to finite temperature!^ It is characterized by i] K = rj' K , 
with 

r) K (d c ,D)=r 1 r :e (4d c ,D), (132) 

where ryP uro (d c , D) is the SCSA exponent for r) K found in Sec 12 . 4 . 2l characterizing a 
homogeneous polymerized membrane at a finite temperature. The underlying rea- 
sons for this amazing connection between roughness exponents at the disorder- and 
thermally-dominated fixed points is unclear. However, this SCSA prediction agrees 
with the l/d c - and e-expansions to lowest order in respective small parameters. For 
a physical membrane, D = 2, d c = 1, SCSA predicts: 

r} = 2/(2 + V6) = 0.449 , C = 0.775 , (133) 

that compares quite well (and much be tter than the lowest order e- or l/d c - 
expansions) with the numerical simulatiorl^ result £ = 0.81 ± 0.03 for a heteroge- 
neous polymerized membrane. 

3.3.2. Long-range disorder 

Above analysis of short-range impurity disorder can be easily extended to treat 
disorder with long-range correlations, that can arise from weakly correlated dis- 
tribution of frozen-in dislocations and disclinations, random grain boundarieJ^, 
and from polymerized-in quasi-long-range correlated lipid tilt (or other membrane 
vector) order. At long scales, such disorder can be characterized by variances with 
power-law Fourier transforms: 

A c (q) = A c q- Z < , (134) 

Ai, a (q) = Ai.ag-* 1 -' , (135) 

where z c and z% t 2 are curvature and stress disorder correlation exponents. Such 
long-range disorder considerably enriches the phase diagram of heterogeneous poly- 
merized membranes, introducing a number of new flat-glass phases, that are sum- 
marized as function of value of these range exponents in Fig^J For sufficiently 
short-ranged disorder (both z c and small), and for finite and zero temperature, 
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we respectively recov er the SCSA exponents for the Aronovitz-LubenskjEEED anc j 
Morse-LubenskjE3Hl fixed points. 
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Fig. 11. Domain of stability of the flat phases as a function of Z/j,,z K . (1) Disorder-dominated 
phases (£' > £): long-range stress glass (LRSG), long-range curvature glass (LRCG), long-range 
mixed glass (LRMG). (2) Temperature-dominated phases (f < f): LR curvature (LRCAL), LR 
mixed (LRM) and LR stress (LRS). (3) LRSM: marginal phase with £ = The shaded area 
corresponds to a region of thermal phase transitions between several stable phases (LRCG and 
the others). The region where the membrane crumples is indicated. 



More generally, the nature of the stable phases strongly depends on the value 
of the disorder-range z exponents and divides into three classes: (1) £ > C' with 
temperature dominated roughness, (2) C < C with disorder dominated roughness, 
and (3) £ = C with equal scaling of the disorder and thermal contributions to the 
membrane roughness. Each one of these regions then further subdivides into distinct 
flat-glass phases depending on whether long-range curvature, stress, or both types of 
disorders are relevant. That is, in the presence of long-range disorder four new flat- 
glass phases, stable at finite temperature appear: (i) short-range (SR) curvature and 
long-range (LR) stress disorder (LRSG in Fig lTT|l . (ii) LR curvature and SR stress 
disorder (LRCG), (hi) LR disorder in both curvature and stress disorder (LRMG), 
and (iv) zero curvature disorder and LR stress disorder (not represented in Fig llljl . 
In addition to these flat-glass phases three corresponding temperature-dominated 
flat phases (LRS, LRCAL, LRM), and two phases for which £ = C appear. In the 
shaded area in Fiell II several of these phases are stable. Phase transition controlled 
by strength of disorder and/or temperature is therefore expected between them. 
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For sufficiently long range disorder correlations (large z's), the dominant rough- 
ness exponent (f reaches 1, presumably indicating disorder-driven crumpling tran- 
sition, and therefore breakdown of the weak-disorder expansion about (on average) 
flat phase. The expected energy-driven crumpling instability to a qualitatively dis- 
tinct isotropic "crumpled-glass" state is a priori of entirely different nature than the 
entropy-driven crumpling transition predicted for phantom membranes. 



3.4. Strong quenched disorder: "crumpled-glass" 

Description of the strongly disordered crumpled-glass phase and the associated 
transition is significantly more complicated, because in addition to complexi- 
ties of conventional spin-glasses, nonlocal self-avoiding interaction must be in- 
cluded. Such crumpled phase is characterized by a vanishing average tangent 
field (d a ri) (hence crumpled), but with a nonzero crumpled-glass order parameter 
(d a ri){dprj) analogous to Edwards- Anderson spin-glass order parameter in disor- 
dered magnets!^ Some progress toward description of such crum pled -glass phase 
in phantom membranes was made by Radzihovsky and Le Doussa^ , by utilizing 
a l/rf-ex P ansion! 41 l 42 l 43 l 

In the limit of large embedding dimension, d — > oo, the homogeneous mem- 
brane model can be solved exactly.^ In contrast, for a disordered membrane even 
in the d — > oo limit the exact solution of the crumpled-glass phase appears to be 
intractable. The difficulty arises from the tensor structure of the crumpled-glass or- 
der parameter, that leads to a problem of matrix field theory, a notoriously difficult 
problem. However, some progress can be made within an additional mean-field like 
approximation that ignores fluctuations in the tensor crumpled-glass order param- 
eter. To simplify technical aspects of the presentation it is convenient to specialize 
to a purely scalar stress-only disorder, described by random Gaussian, zero-mean 
dilations and compressions in the locally preferred metric, </L(x) = S a f3(l + Sg°(x)). 
A much more questionable (but technically necessary) approximation is omission of 
the self-avoiding interaction, that is undoubtedly important in the crumpled-glass 
phase. 

Because of the isotropic nature of the crumpled-glass phase, the model must 
be formulated in terms of d-dimensional conformation vector r(x). The effective 
Hamiltonian is given by 



F[r\=d d u x 



f |V 2 rf + £{d a ? ■ d p r- ^(x)) 2 + ±(d a ? ■ d a r- gU^Y 



(136) 



where the elastic moduli were rescaled by d so as to obtain sensible and nontrivial re- 
sults in the limit d — ► oo. Replicating F allows averaging over quenched disorder!^ 
Then, introducing two Hubbard-Stratanovich fields Xa/3 and Qabafi ^° decouple 
replica diagonal and off-diagonal nonlinearities, respectively, leads to an effective 
Hamiltonian that is quadratic in f. This allows formal integration over r, that is 
conveniently done around background configuration rp. Now, ignoring fluctuations 
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in the Hubbard-Stratanovich fields, the values of the order parameters d a fo and 
QabafSij are determined by minimizing the resulting replicated free energy, together 
with the equation of constraint relating Xaap to these order parameters. 

Assuming that the replica symmetry breaking does not o ccur u ntil higher order 
in 1/d, as it happens in the random anisotropy axis 

mode J57|58| 

we look for the 

saddle point replica-symmetric solution of the following form, 

f% = (x a e a , xZ a p = X s <*p , Qlbapij ■= Q8 a pSij(l - S ab ) . (137) 
The corresponding saddle-point equations for £, \ an d Q are given by: 

2 AQ/T 




T 


f A d D k 


2D 


Jo i^) D 






Q ^ 






0. 





± 2 + X+£Q (k^+x+^Q) 2 



A ,D 



d u k 1 



o (2tt) b ( K p + x + 4 Q) 2 



,(138) 
(139) 
(140) 



where A is scalar stress disorder variance and A = (2/x + .DA)A. For the special 
case of homogeneous membranes, A = 0, these equations reassuringly reduce to 
those found by David, et al.El and describe the crumpled-to-flat transition exactly 
to leading order in 1/d. 

For a heterogeneous membrane (A > 0) there are three distinct solutions to 
these saddle-point equations, corresponding to three different possibilities for the 
values of the pair of order parameters £ and Q. 

C = 0, Q = , (141) 
C ? , Q^0 , (142) 
C - , Q^0 . (143) 

that correspond to the crumpled phase, flat phase and crumpled-glass phase of the 
membrane, respectively. 

Critical properties of these three phases and phase bou ndar ies between them can 
be obtained from a straightforward analysis of Eas ll38ll4()[ ^lThe phase behavior is 
summarized in Fig ll2l illustrating that within this approximation the lower-critical 
dimension for the flat phase in the presence of quenched disorder is Df^ = 4, and 
therefore for this model (in d — > oo limit) only crumpled-glass and thermal crumpled 
phases survive in D — 2 membranes. 

For the flat phase saddle point equations give 

f-^-S) ('-£)• • (144) 

where critical crumpling transition temperature T c and critical value of disorder A c 
(not to be confused with the variance of the curvature disorder from Secs l3. 213.3. 21 . 
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defining the rectangular boundaries of the flat phase, are given by 

i_ r A _d^±_ a = j_ r A rf D fc 1 

£> Jo (2tt) d nk 2 ' c 2D 



(27T) 



(145) 



and A 1 = 1 + (a + Df3)A/T. The power-law vanishing of these order parameters 
according C ~ (T c - T)4, ( ~ (A c - A)' 3 ! and Q ~ (T c - T)' 3 ?, defines the 
corresponding /? exponents: = 1/2, /j£ = 1/2, /3® = 1. 

Outside this rectangular region £ vanishes and the membrane undergoes a crum- 
pling transition out of the flat phase. For A < A c and as T — > T c the transition is 
to the crumpled phase, while for T < T c and as A — > A c flat phase is unstable to 
the crumpled-glass phase. 



D>4 

crumpled glass, 
Q*0, C=0 




crumpled 

Q=o, C=o 



o 



D<4 

crumpled glass 
Q*0, C=0 

'crumpled 

Q=0, C=0 



flat t c 

Q*=0, <j0 (b) 



Fig. 12. A phase diagram for a disordered polymerized membrane showing (a) crumpled, flat- 
and crumpled-glass phases for D > 4. (b) In the physical case D = 2 < 4 the flat phase does not 
appear in this d — + oo limit, but is expected to when 1/d corrections are taken into account. 



The glass susceptibility near the transition from the flat phase to the crumpled- 
glass phase can also be easily calculated by introducing an external field hij a p = 
hSijS a /3 conjugate to the crumpled-glass order parameter Qabapij- The result- 
ing saddle-point equations then lead to the crumpled-glass susceptibility, \ sg — 
dQ/dh ~ (A c - A) -7.*" , with lsg2 = 1 . 

Similarly, upon approach to the flat phase from the crumpled-glass (character- 
ized by Q ^ and ( = 0) the tangent susceptibility X( diverges as tangent order ( 
spontaneously develops. Turning on an external field / that couples to the tangent 
order parameter, leads to X( =d(/df~ (A — A c )~ 7 < 2 giving 7^2 = 2/|4 — D\, that, 
as expected diverges at the lower-critical dimension Dfi. = 4 of the flat phase. 

Finally, we look at the transition between the crumpled and crumpled-glass 
phases. The crumpled-glass susceptibility \sg near this transition is given by \sg ~ 
(A C (T) — A)~ 7ssl , with 7„ 9 i = 1, as at the flat-to-crumpled-glass transition, except 
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for the modified phase boundary that is nonzero for any D: 



A-\T) = 



d D k 



2DJ (2n) D («jfc2+ x (T)r 



(146) 



and together with saddle-point equations and Q = 0, also defines the phase bound- 
ary between the crumpled and crumpled-glass phases for D > 2, A C (T) — A c ~ 
(T - T c )*, with 4> the crossover exponent <f> = \D - 4|/(D - 2). 

As noted above, d — ► oo analysis predicts an instability of the flat phase of 
D = 2 membranes to any amount of disorder (D^. = 4). A computation of 1/d 
corrections for the disordered membrane is technically quite challenging and remains 
an open problem. However, quite generally, anomalous elasticity generated by 1/d 
corrections (e.g., finite rj K = 0(1/ d) exponenlP) strongly suggests the lowering of 
Dfi: down to Dy. = 4 — 0(l/d). This is supported by the Harris criterion applied to 
the buckling transition, controlled by the Aronovitz-Lubensky fixed point, that leads 
to stability of the flat phase (and the AL fixed point) as long as r] u remains positive. 
This is also consistent with the e = 4 — £>-expansion analysis of Radzihovsky and 
Nelson (performed for arbitrary (fpHI discussed in Sec 13. 3. il that the lower-critical 
dimension is reduced down to Dfc = 4 — 4/d, A phase diagram for D < 4 consistent 
with the nature of the 1/d corrections is illustrated in FigllHI 







D<4 

crumpled glass, 
Q*0, ^=0 




Fig. 13. Conjectured phase diagram for a disordered membrane D < 4 when 1/d corrections are 
taken into account, that allow a region of flat phase of size 1/d to appear. 



Finally we observe that the crumpled-glass phase can be destroyed by apply- 
ing an external tension to the membrane's boundaries. The metastable degenerate 
ground states would disappear and the average of the local tangents would no longer 
vanish. In this respect an external stress would be analogous to an external magnetic 
field in spin systems. As the stress is reduced the membrane would slowly return to 
the glassy phase but with some hysteresis. The line separating the regions of stable 
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and metastable degenerate states is then the analogue of the d'Almeida-Thouless 
studied in great detail for the real spin-glasses!"-^ 



4. Interplay of Anisotropy and Heterogeneity: Nematic Elastomer 
Membranes 

We would like to conclude these lectures with a discussion of a new exciting realiza- 
tion of polymerized membranes, nematic elastomer membranes .UJ The motivation 
for their study is driven by r ecen t experimental progress in the synthesis of ne- 
matic liquid-crystal elastomer^2l statistically isotropic and homogeneous gels of 
crosslinked polymers (rubber), with main- or side-chain mesogens, that can sponta- 
neously develop nematic orientational order, accompanied by a spontaneous uniaxial 
distortion illustrated in FiglT4l 





Nematic 



IN 



Isotropic T 
Q=0, u=0 



Fig. 14. Spontaneous uniaxial distortion of nematic elastomer driven by isotropic-nematic tran- 
sition. 



Even in the absence of fluctuations, bulk nematic elastomers were 
and later observed to display an array of fascinating phcnomcn 



redictecPSl 

the most 

striking of which is the vanishing of stress for a range of strain, applied transversely 
to the spontaneous nematic direction. This striking softness is generic, stemming 
from the spontaneous orientation al sy mmetry breaking by the nematic stat e I I, 
accompanied by a Goldstone modj^j that leads to the observed soft distortion and 
strain-induced director reorientation^, illustrated in Fig^J The hidden rotational 
symmetry also guarantees a vanishing of one of the five elastic constanta^S that 
usually characterize harmonic deformations of a three-dimensional uniaxial solidP^I 
Given the discussion in Sec l2.4.ll not surprisingly, the resulting elastic softness leads 
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to qualitative importance of thermal fluctuations and local het erogeneity. Similar to 
their effects in smectic and columnar liquid crystals 05211011121 thermal fluctuations 
lead to anomalous elasticity (universally length-scale de penden t elastic moduli) in 
bulk homogeneous elastomers with dimensions below 3ISSEI anc [ below 5, when 
effects of the random network heterogeneity are taken into account!^ 




(a) (b) 

Fig. 15. (a) Simultaneous reorientation of the nematic director and of the uniaxial distortion is a 
low-energy nemato-elastic Goldstone mode of an ideal elastomer, that is responsible for its softness 
and (b) its flat (vanishing stress) stress-strain curve for a range of strains. 



This rich behavior of bulk elastomers provided strong motivation to study ne- 
matic elastomer mem branes (D-dimensional sheets of nematic elastomer fluctuating 
in d dimensions) A model of such a membrane must incorporate both network 
anisotropy and heterogeneity discussed in previous sections. However, an important 
distinction from explicitly anisotropic membranes discussed in Sec|21is that the ne- 
matic anisotropy is spontaneously chosen in the amorphous (initially statistically 
isotropic) elastomer matrix. At harmonic level this in-plane rotational symmetry 
can be captured by a two-dimensional harmonic effective Hamiltonian 

K° NE = \j d 2 x[ Kxx (d 2 h) 2 + Kyyitfh) 2 + 2 Kxy (d 2 x h)(d 2 h) + K y (d 2 u x ) 2 + K x (d 2 u y ) 2 

~^~^x ij^xx ) ^yi^-yy*) ^^xy^xx^yyj j (1^^) 

with e aa — d a u a (no sum over a implied here), and characterized by a uniaxial 
phonon elasticity with a vanishing transverse shear modulus, fi xy . This latter feature 
is what distinguishes a nematic elastomer membrane from an explicitly anisotropic 
membrane discussed in Sec|21 The vanishing in-plane shear modulus captures at the 
harmonic level the invariance of the free energy under infinitesimal rotation of the 
nematic axis and the accompanying uniaxial distortion of the elastomer matrix. To 
ensure an in-plane stability curvature phonon elastic energies are included in 7i5v S . 
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As a result, in the putative flat nematic phase of an elastomer membrane, the 
phonons have qualitatively "softer" harmonic elasticity than in a conventional poly- 
merized membrane discussed in SecEl Consequently, as in other "soft" systems 
(e.g., smectic and columnar liquid crystals phases; see discussion in Sec l2.4.f|) . in 
the presence of thermal fluctuations and/or heterogeneities, nonlinear elastic terms 
are essential for the correct description. 

First principles derivation of the nematic elastomer model, that incorporates 
(hidden) in-plane rotational invariance at nonlinear level is somewh at invo lved and 
we refer an interested reader to recent detailed work on this subject E *6 | 37 | j n snor ^ 
one starts out with a model of a statistically homogeneous and isotropic elastic mem- 
brane coupled to a nematic in-plane order parameter Q a p. The rotational symmetry 
is then spontaneously broken by the nematic ordering at the isotropic-nematic tran- 
sition, that also induces a spontaneous uniaxial in-plane distortion of the elastomer 
matrix, characterized by a strain tensor vP a n. Expansion about this flat uniaxial 
state, ensuring underlying in-plane and embedding space rotational invariance leads 
to a nonlinear clastic Hamiltonian of a nematic elastomer membrane. Its form is 
that of the 'H% E) Eg 11471 but with the harmonic strain e a p replaced by a nonlinear 
strain tensor w a p, that incorporates both in-plane and height nonlinearities of the 
form (d x u y ) 2 and (d x h) 2 , respectively. 

In D = 2 both phonon and height anharmonic terms are strongly relevant (their 
perturbative corrections grow with length scale) and therefore must be both taken 
into account. One approach is to generalize the above model to a .D-dimensional 
nematic elastomer membrane and perform an RG calculation controlled by an e- 
expansion. However, it is easy to show that upper-critical dimensions for these non- 
linearities are different, with height undulations becoming relevant below D\ c — 4, 
and smectic-like and columnar-like in-plane nonlinearities with upper-critical di- 
mensions of — 3 and D^J = 5/2, respectively. Hence for D > 3, in-plane 
phonon nonlinearities can be neglected, with height undulation nonlinearities (of 
the type studied in Secs l2.4.2l2.5.1l the only remaining relevant ones for D < 4, and 
controllable close to D = 4 with an e = 4 — _D-expansion. 

Generically, one would expect these undulation nonlinearities to renormalize 
bending rigidities well as in-plane elastic moduli \ a p, leading to anomalous 

elasticity with K a p(q) ~ <7~'' K , X a fi(q) ~ q Vx . As discussed in Sec l2.4.2l here too, 
rotational invariance imposes an exact Ward identity between exponents: 

2r} K + T) X = e = 4 - D. (148) 

However, it is not difficult to show^J that once in-plane nonlinearities are neglected 
(legitimate for D > 3), the harmonic phonons can be integrated out exactly, and 
lead to a purely harmonic effective Hamiltonian TL[h\. Therefore there is a strict 
nonrenormalization of the bending rigidity tensor n a p for D > 3. This together 
with the Ward identity, Ea 11481 gives 



ri K = 0, rix = 4 - D , 



(149) 
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a result that is supported by a detailed renormalization group calculation!^ This 
analysis also makes contact and recovers some of the results previously obtained 
in the studies of isotropic polymerized membranes. In particular, the previously 
seemingly unphysical, the so-called "connected fluid'^3 is realized as a fixed point of 
a nematically-ordered elastomer membrane that is unstable to the globally stable 
nematic-elastomer fixed point 

Despite of some success, there are obvious limitations of above description, most 
notably in its application to the physical case of D = 2 elastomer membranes and 
inclusion of the (usually more dominant) local network heterogeneity. The first 
shortcoming primarily has to do with the neglect of in-plane elastic nonlinearities, 
which, near the Gaussian fixed point become relevant for D < 3. While it is very 
likely that the subdominance of these in-plane nonlinearities relative to the undula- 
tion ones will persist some amount below D = 3 1^21 we expect that in the physical 
case of D — 2 all three nonlinearities need to be treated on equal footing. Carrying 
this out in a consistent treatment remains an open and challenging problem. 

Secondly, elastomers are only statistically homogeneous and isotropic, exhibit- 
ing significant local heterogeneity in the polymer network. As we saw in Sec|31 such 
internal quenched disorder has rich qualitative effects in ordinary polymerized mem- 
branes. Furthermore, recent work by Xing and Radzihovsky has demonstrated^^ 
that interplay between nonlinear elasticity and random strains and torques (due 
to network heterogeneity) leads to disorder controlled anomalous elasticity even in 
three-dimensional bulk nematic elastomers. Because nematic elastomer membranes 
are far softer than ordinary polymerized membranes and their bulk analogs, we 
expect network heterogeneity to have strong and rich effects in these systems. Con- 
siderable research remains to work out the resulting phenomenology. 




1 FT 



1 1N 



l CF 



Fig. 16. A possible phase diagram for ideal nematic elastomer membranes. As the temperature is 
lowered a crumpled membrane undergoes a transition to isotropic flat phase at Tq p , followed by a 
2D in-plane isotropic-nematic like transition to an anisotropic (nematic) flat phase. As T is lowered 
further, this anisotropic flat phase becomes unstable to a nematic tubule phase, where it continu- 
ously crumples in one direction but remains extended in the other. At even lower temperature, a 
tubule-flat transition takes place at Ttf- 



It is interesting to conclude with a general discussion of the global conforma- 
tional phase behavior of nematic elastomer membranes. As with ordinary polymer- 
ized membranes, upon cooling, isotropic elastomer membranes should undergo a 
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crumpling (flattening) transition from the crumpled to flat-isotropic phase. Upon 
further cooling, an in-plane (flat) isotropic-to-(flat) nematic transition can take 
place, leading to a flat membrane with a spontaneous in-plane nematic order. As 
for explicitly anisotropic membranes discussed in Sec[5J such nematically-ordered 
elastomer membranes should undergo further transition to a nematically-ordered 
tubule phase. Because of the in-plane rotation symmetry that is spontaneously (as 
opposed to explicitly) broken, such nematic tubule will exhibit in-plane elasticity 
("soft" phonons) that is qualitatively distinct from tubules discussed in Sec l2.5.ll 
and will constitute a distinct phase of elastic membranes!^ One interesting sce- 
nario of phase progression is the nematic-flat to nematic-tubule to nematic-flat 
reentrant phase transitions, driven by competition between growth of nematic order 
(anisotropy) and suppression of membrane's out-of-plane undulations upon cooling, 
as schematically illustrated in Fig^5] Considerable research to elucidate the nature 
of the resulting fascinating phases and transitions remains .^31 

5. Summary 

In these notes, I have presented a small cross-section of the beauty and richness 
of fluctuating polymerized membranes. I have demonstrated the importance of the 
in-plane order in determining the long-scale conformations of these clastic sheets, 
by discussing in-plane anisotropy and local random heterogeneity and showed that 
these lead to a rich and highly nontrivial phenomenology. 
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